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Abstract I briefly review the method of population 
synthesis of binary stars and discuss the preliminary 
results of a study of the Galactic population of subd- 
warf B stars. In particular I focus on the formation 
of (apparently) single sdB stars and their relation to 
(apparently) single helium-core white dwarfs. I discuss 
the merits of mergers of two helium white dwarfs and 
interactions with sub-stellar companions for explaining 
these single objects. A preliminary conclusion is that 
the current observations suggest both mechanisms may 
contribute, but that the helium white dwarfs are likely 
formed in majority from interactions with sub-stellar 
companions. 



1 Introduction 

In a volume like this one, no lengthy introduction to 
subdwarf B (sdB) stars is needed, as many aspects 
of their nature are discussed in the different contribu- 
tions. Instead I will briefly describe the motivation to 
study the Galactic population of sdB stars. Because 
rather special circumstances are needed to form an sdB 
star (the hydrogen envelope needs to be expelled at the 
right moment), they offer one of the best populations 
to constrain binary-evolutio n models, in particular the 
common-envelope phase (e.g. Han et al. 2003t Hu et al. 
200i . 



Secondly, the fact that some sdB stars are pul- 
sators allows asteroseismological determination of ac- 
curate masses and internal structure, which provides 
tests of stellar evolution and in particu lar an excellent 



way to constrain their formation (e.g. iFontaine et 



boost and contributions of Van Grootel and Hu in this 
volume) . 
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Finally, the relative brightness of sdB stars allows 
us to study them to much larger distances than the 
other objects that probe the late phases of evolu- 
tion of low- and intermediate-mass stars, the white 
dwarfs. Indeed, they are likely responsible for much 
of the UV light observed in elliptical ga laxies (e.g. 
Han. Podsiadlowski. and Lvnas-Gravl 120071 and Han, 
this volume). 

However, I think the close link of sdB stars with 
white dwarfs - either low-mass helium-core white 
dwarfs as "failed" sdB stars, or the low- mass C/0 core 
white dwarf descendants of sdB stars - deserves more 
attention. In the remainder of this paper I will discuss 
the principle of population synthesis (Sect. [2]), a model 
for the Galactic population of sdB stars (Sect. ^ a dis- 
cussion of sdB stars in relation to white dwarfs in Sect. 2] 
and the probability of finding sdB stars in binaries with 
neutron stars or black holes (Sect. [5]). 



2 Population synthesis 

One of the ways to constrain the binary evolution us- 
ing sdB stars, is using so called population synthe- 
sis (for a general descrip tion of th e different popula- 
tion synthesis codes see iKornilov a nd Lipunov 1983; 
Dewey and Cordea 1987 : Lipunov aiid Postnovl Jl988 : 
de Kool"l99d;'Yungelson and Tutukov"l99l';'Kolb''l993l 

Pols and Marinus 1994- Portegies Zwart and Verbunti 

1996 : Vanbeveren. van B ever and de Dondeilll997 : Frver. Burrows, anc 



1998: 



Nelemans et al\ 



2001 



Hurlev. Tout, and Pols 



2002t iHan et al\ 120031 : iBelczvnski et al\ l2008f) . The 
basic concept of this method is that a large number 
of binary evolution scenarios are calculated using (ap- 
proximations to) the evolution of stars (both single and 
those affected by binary interactions), in combination 
with modelling or recipes of binary interactions. The 
scenarios, or their relative weighting are chosen in such 
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a way as to represent assumed initial parameter distri- 
butions (initial primary and secondary mass and orbital 
period). The outcome of this procedure yields a model 
for the Galactic population of binaries, that only needs 
to be normalized in a suitable way to the estimated 
total number of (binary) stars in the Galaxy. 

When studying objects using population synthesis 
one should keep in mind that several crucial ingredi- 
ents of the models are poorly known. In particular 
the effect of binary interaction and mass transfer on 
the further evolution of the system, the strength of the 
mass loss (in particular for massive stars), the effect 
of supernova explosions, and evolution of peculiar stars 
that can only be formed via binary interactions. In ad- 
dition, in order to quantitatively compare the models 
to observed populations, the initial parameter distribu- 
tions must be known. The distribution of masses for 
the most massive (primary) component of the binary 
is usually taken to be the initial mass function (assum- 
ing the effect of binaries on the determination of the 
IMF is dealt with already). For the secondary mass, 
the observational constraints are of course heavily bi- 
ased towards similar masses, as otherwise the light of 
the lower-mass companion (secondary) is buried by the 
bright primary. Finally the distribution of orbital pe- 
riods must be patched together from spectroscopic and 
eclipsing binaries for the shorter periods, via visual bi- 
naries to common proper motion pairs for the wider 
syst ems. Gurrently most people use a standard IMF 




Fig. 1 Representation of our assumed Galactic s tar- 
formation model (based on lBoissier and Prantzoj|l999l ) for 
three different Galactic ages. The star formation proceeds 
from the center, including the Bulge, to the outer disk. The 
thick disk and Halo are not modelled. 



Fig. [T] we show three snap shots of the formation of the 
Galaxy, including Bulge and thin disk. The thick disk 
and Halo are not included in the model. 



3 Galactic populations of sdB stars 



a flat mass ra- 


Han et al. 


( 


2002. 


2003 


) and 


Yungelson and Tutukov 



tio distribution (i.e. for a given primary mass, each 
mass lower has equal probability) and a initial semi- 
major axis distribution that has equal numbers of sys- 
tems in logarithmic intervals from very close binaries to 
typically ^10^ AU. The fact that most people use the 
same initial conditions is good for comparing results, 
but the currently used distributions are most likely too 
simple. 

When considering low- and intermediate-mass stars 
also the normalization is an issue to be considered. In- 
teresting phenomena or types of systems (such as the 
formation of sdB stars) may occur several Gyr after 
the formation of the main sequence binary. There- 
fore the star formation history of the Galaxy be- 
comes important. In most studies the star forma- 
tion has been taken to be constant or, for old pop- 
ulations, to be a single star burst long ago. In 



sis of sdB stars, building on many earlier studies of 
the forn iation of helium core burning stars in bina - 
ries (e.g. lWebbin"^ll984[[Tutukov and YungelsoJl990l ). 



So why is it useful to do another population synthesis 
study? Firstly, because of all the uncertainties in the 
binary evolution and assumptions about initial param- 
eters and star formation, it is good to compare simi- 
lar calculations from different groups. In particular, we 
have concluded from our investigations of the formation 
of double white dwarfs that a different description for 
the common-envelope phase was needed in our model 
j Nelemans et al. 2000h. a lthoug h this has been used by 
Yungelson and Tutukovl (j2005h as well. Secondly, as 



Nelemans. Yungelson. and Portegies ZwartI (|2004[ ) we 
started to work with a more realistic star-formation 
history, based on the gal axy formation models of 
Boissier and Prantzos This is an inside-out star 

formation model, that peaked early in the history of the 
Galaxy and has decreased substantially since then. In 



discussed above, we use a rather different star forma- 
tion history. 

The model p resented here is basically the same as 
that d iscussed in Nelemans. Yungelson. and Portegies Zwart^ 



1 2004^ . with the exception that for our alternative 
common envelope we use a value of 7 = 1.5, rather 
than 1. 75. This value fits most of the observed bi- 
naries ([Nelemans and ToutI 120051 ) . Note that we still 
apply the alternative common envelope only to the 
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cases where none of the two stars is a compact ob- 
ject and the common envelope is triggered by dy- 
namical unstable Roche-lobe overflow (rather than a 
tidal instability). The rea son is that, even though in 
Nelemans and ToutI (|2005l ) we conclude that the al- 
ternative formalism seems to be able to explain all 
observed binaries, the motivation for the alternative 
formalism is the large amount of angular momentum 
available in binaries with similar mass objects. 

The basic stellar evolutio n inpu t is described in 



Table 1 Local Space densities and total numbers in the 
Galaxy of different classes of sdB stars and binaries 



Portegies Zwart and Verbunti (|1996l ): iNelemans et al. 
( 20011 ) . The stellar evolution tracks are somewhat out- 
dated and we are in the process of updating them. 
For the moment the evolution of the core of the star 
is determined by core-mass luminosity relations and 
following of the growth of the core by shell burning 



( Nelemans et al\ 120011 ). In our population synthesis 
calculations, we classify all helium core burning stars 
that have lost their hydrogen envelope as sdB stars. 

The only update for this work is that we allow degen- 
erate helium cores that get exposed via binary interac- 
tions and have a core mass close to the mass where the 
heliiini flash happens (within '^0.02 Mq, see lHan et al. 
200i IHu et oil I2OO7I) to ignite and become sdB stars 
(|D'Cruz et aLlll996l) . This feature was not included in 
our models before, as those cores were assumed to be- 
come helium white dwarfs. Because the helium flash for 
single stars in the stellar models that we use happens at 
a rather low core mass of 0.446, this channels produces 
sdB stars with masses between 0.426 and O.446M0. 

3.1 Results 

I present here the preliminary results of a model for 
the Galactic population of sdB stars that will be pub- 
lished in more detail elsewhere. We distinguish three 
classes of objects, sdB stars with main sequence (MS) 
companions, with white dwarf (WD) companions and 
single sdB stars. 

The number of systems in the Galaxy at present 
is given in local space densities in Table [U which is 
determined by the time and position dependent star- 
formation rate. We give several different values for dif- 
ferent selection effects that may play a role when com- 
paring with observations. First for sdB stars with MS 
companions, we select those where the sdB star is at 
least as bright as the companion in the V-band. Sec- 
ondly, for both WD and MS companions, we give the 
density of systems with periods shorter than 10 days, 
which is about the longest period systems that can be 
found relatively easily using spectroscopy (e.g. Geier, 
this volume). 

The total space density of sdB -|- MS binaries is 
about a factor of 5 larger than that of sdB -I- WD 
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Fig. 2 Period - main sequence mass distribution for the 
sdB -I- MS binaries with V-band magnitude of the binary 
brighter than 20 but no further selection effects. The dif- 
ferent formation channels producing this distribution are 
discussed in the text. 



binaries, but for the "observable" systems, the space 
densities are comparable. The space density of double 
helium white dwarf mergers is roughly the same as the 
sum of the the "observable" binaries. For comparison 
with other population synthesis studies, I show in the 
last column of Table [T] the total number of systems in 
the Galaxy. 

In Fig.[2]I show the distribution of periods and com- 
panion masses for the sdB -f- MS binaries (without the 
selection effects mentioned above, but for a magnitude 
limited sample). This gives a n ice view on the differ- 
ent formation channels (see also lHan et a/.l 120021 ). The 
patch at the top, at relatively short periods and mas- 
sive companions is a combination of common-envelope 
ejection of stars with non-degenerate cores, where rel- 
atively large separations become a lot shorter (darkest 
region at shortest periods) , and stable Roche-lobe over- 
flow, which in our model can only happen in relatively 
short period binaries, when the giant fills its Roche lobe 
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Fig. 3 Mass distribution of sdB stars in close binaries. 
Dark liistogram (total): sdB + MS stars (with KidB < Vms)- 
The solid histogram gives the subset of systems with low- 
mass (<0.15 Mq) companions. Light histogram: sdB + WD 
systems. 



in the Hcrtzsprung gap, widening the orbit (wider pe- 
riods in that patch). However, all these systems are 
"unobservable" as the condition T^dB < Vms effectively 
limits the mass of the MS companion to below ~0.8 
Mq. The patch at the right, at very long periods, is 
due to our alternative common-envelope prescription 
for stars near the tip of the red giant branch, leading 
to mild or hardly any spiral-in. Finally, the triangular 
patch at the bottom is for stars that enter a common 
envelope via tidal instability, which we treat with the 
energy balance in our model (see also Soker, this vol- 
ume) . The low- mass and short period part of that patch 
could be associated with the HW Vir type binaries of 
which th ere are sever al observed with periods around 
0.1 d (see lHebeillioogI) . The periods on our model are a 
bit wider (and less concentrated) than those observed, 
which may indicate that we assume a too large effi- 
ciency in the common envelope (we assume aX = 2). 
A detailed comparison of this model with the observed 
systems will given in the forthcoming publication. 

In Figure |3] I show the mass distribution of close bi- 
nary (periods less than 10 days) sdB stars with MS 
(dark) and WD (light) companions. For the systems 
with MS companions I only show systems where the es- 
timated V-band magnitude of the sdB stars is smaller 
or equal to that of the MS star. The dark solid his- 
togram gives the subset of systems with very low-mass 
(<0.15 Mq) companions, to be compared to the num- 
ber of HW Vir stars. 

The sdB stars with MS companions all come from bi- 
nary interactions close to the tip of the red giant branch 
and thus have masses between 0.426 and O.446M0. The 



Fig. 4 Mass distribution of the merger products of double 
helium white dwarfs. The solid histogram below 0.4 Mq 
shows those merger products that to not have enough total 
mass to ignite helium and become an sdB star. Above 0.4 
M0 the ones that do ignite and thud become single sdB 
stars are shown (open histogram). 



ones with WD companions can also originate in more 
massive stars (M > 2Mq) that ignite hehum non- 
degenerately in their cores and thus can both be lower 
as well a higher mass than those from the tip of the 
giant branch. Observable sdB stars with MS compan- 
ions cannot originate from stars with non-degenerate 
cores in our models, because these stars are typically 
more compact and, because of their larger total mass, 
have more massive envelopes. They could only be un- 
bound by massive MS companions that would outshine 
the sdB star. In Fig. 2] the distribution of single sdB 
stars that originate from the merger of two helium white 
dwarfs (see Sect. |4]). The result of mergers are, of 
course, in general more massive. 



4 Single sdB stars and single helium-core 
white dwarfs 

Here I further discuss the formation of (apparently) sin- 
gle sdB/sdO stars. Different models have been pro- 
posed: enhanced m ass loss on the red giant branch 
( D'Cruz et aZlll996h. t he in erger of two helium white 



dwarfs (e.g. Webbinlj Il984l for sdO stars) or inter- 
actions with low-mass (sub-stellar) companions (.Sokeil 
1998h . In the remainder we discuss the latter two, as it 
is unclear if the enhanced wind scenario actually hap- 
pens. 

Interestingly, in both these cases the same mecha- 
nism will also produce helium white dwarfs. In the case 
of the merger when the total mass is too low to ignite 
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Fig. 5 Outcome of the interaction of a 1 M0 giant witii a 
sub-stellar companion as function of the initial period and 
mass of the companion (top). Above the solid line the sub- 
stellar companion survives, below it evaporates inside the 
giant envelope. Between the dotted and solid line a signifi- 
cant fraction of the envelope mass is lost and the rest may be 
lost by stellar winds. The lower panel gives the mass of the 
core at the moment of the interaction, i.e. the mass of the 
resulting stars. In order to get an sdB star the core needs 
to be within a very small mass region (see Sect.[3)l, so only 
systems with initial periods between 1200 and 1300 days can 
produce sdB s tars, the rest produces a hel ium white dwarf. 
Adapted from iNelemans and Tauri^ (| 19981 ) 



the h elium (probably around 0.4 Mq, see iHan et al. 
2002!) ■ and in the case of an interaction on the giant 
branch when the mass of the degenerate core is too low 
to ignite helium after the hydrogen envelope is lost (as 
proposed for the observed s ingle helium white dwarfs in 
Nelemans and Taurij 19981) . In both cases the ratio of 
the formation of single helium white dwarfs and single 
sdB stars may be different, and these model ratios can 
then be compared to the observed ratio of single helium 
white dwarfs to single sdB stars. 

In Fig. |4] I show the mass distribution of the dou- 
ble helium white dwarf merger products which shows 
that in the merger case the vast majority ('^90%) of 
objects would ignite helium and thus become sdB (or 
sdO) stars. Fig. [SJ adapted from INelemans and Tauris 
( 19981 ). shows the outcome of a common-envelope of a 
sub-stellar companion with a 1 Mq red giant for dif- 
ferent initial orbital periods and companion masses. If 
I look at systems in which the sub-stellar companion 
survives the common envelope and the degenerate he- 



lium core is massive enough to ignite and produce an 
sdB star (denoted sdB+planet in the figure), the ini- 
tial orbital period of the sub-stellar companion is re- 
stricted to a very narrow range, say between 1200 and 
1300 days. Therefore if we assume as initial period are 
equally likely, the probability of forming a helium white 
dwarf is a lot larger than of forming an sdB star. A 
complicating factor is that even if the sub-stellar com- 
panion merges with the giant (region denoted normal 
single WD in the figure), it may cause significant mass 
loss and spin-up the giant enough that later it will lose 
its envelope in a stellar wind, exposing the core either 
when it is still on the red giant branch (yielding a he- 
lium white dwarf or an sdB star depending on the mass 
of the core), or after it contracted to th e horizontal 
branch following the helium flash (see also|SokerJ|l998; 



Politano et aLll2008l ). We therefore included sub-stellar 



companions in our population synthesis calculation by 
simply extending the flat mass ratio distribution to in- 
clude sub-stellar companions, and looked at the ratio of 
the formation of single helium white dwarfs versus sin- 
gle sdB for interactions on the red giant branch. As the 
mass ratio distribution at these masses and the fraction 
of stars with substellar companions are very uncertain, 
we do not use this calculation to determine the absolute 
number of systems, but only to assess the ratio of he- 
lium white dwarf to sdB stars formation. As expected, 
for companions that survive the common envelope, we 
find a ratio of about 10 for the number of helium white 
dwarfs versus sdB stars, but even for the mergers that 
subsequently lose their envelope, we find more helium 
white dwarfs than sdB stars (by about a factor 3). The 
total ratio comes out at 5.5. 

How does this compare with the observational con- 
straints? From the SPY survey ([Napiwotzki et al 



20011 ) we find that 15 out of 636 white dwarfs observed 



are single helium white dwarfs (see Fig. [6]). Assum- 
ing a birth rate of white dwarfs (or pla netary nebu- 
lae) of 3 X 10-^2pc-3yi.-i (jPottaschlll99i ). this results 
in a formation rate of single helium white dwarfs of 
7.5 X 10~^'*pc^'^yr^^. For the single sdB stars we use 
the estimated birth rate of 4 x 10~^^pc~'^yr~^ ( Heber 
1986, ). which together with a binary fraction of about 
50% gives a birth rate of 2 x lO^^^'pc^'^yr^^. Because 
single helium white dwarfs are more common than sin- 
gle sdB stars, they cannot be predominantly formed 
from mergers. Instead it suggests the the single he- 
lium white dwarfs are formed from interaction on the 
red giant branch. If all are formed in that way, the 
same mechanism will produce single sdB stars at a rate 
comparable with the observed rate and thus may be re- 
sponsible for the production of most of the single sdB 
stars. 
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Fig. 6 Effective temperature - surface gravity plot for all 
white dwarfs observed in the SPY survey ( Napiwotzki et ali 
I2OOII ). The squares show the close binaries, the dots the 
single white dwarfs. The solid lines give colling tracks for 
white dwarfs with mass 0.41, 0.29 0.24 and 0.20 Mq from 
lAlthaus. Serenelli. and Benvenutol (|200ll ) 



However, care should be taken with this reasoning, 
as there are more ways to produce single sdB stars or 
single helium white dwarfs which could influence the 
ratio estimated above. In particular, in Sect. 13. il l have 
shown that sdB stars with low-mass companions in wide 
binaries can easily be more than half of the population 
(Table [T]) even for systems where the sdB star domi- 
nates the (blue) light. Current observational radial ve- 
locity studies are not sensitive to these long periods, i.e. 
would not detect the sdB stars as binaries and thus clas- 
sify them as single. In addition, single sdB stars may 
be fo rmed in so-called AM CVn systems fsee lNelemans 
2OO5I for a review) . Accretion of helium onto a helium 



white dwarf can lead to ignition at the edge and an 
outside in helium burning front, very likely turning the 
accreting white dwarf into a helium burning star (but 
probably without H atmosphere). A de tailed study of 
this sc enario is in preparation. Finally, Justham et all 
( 2009f ) have proposed the companions of stars that ex- 
plode as type la supernova to be observable single he- 
lium white dwarfs possibly affecting that side of the 
single helium white dwarf to sdB star ratio. However, 
the birth rate derived for single helium white dwarfs 
above, is about a factor ten above the rate of type la 
supernova explosions from single degenerate scenarios, 
so this scenario is not important for the overall number 
of single helium white dwarfs. 

I conclude that interaction of sub-stellar companions 
may be a significant contributor to the formation of 
single sdB stars, but a full detailed study, including all 
scenarios and proper selection effects is needed to come 
to firm conclusions. 



5 Massive companions? 

One of the most intriguing results of the last years 
is the discovery of apparently massive, compact com- 
panions to sdB stars (e.g. Geier et al., this volume 
and references therein). These masses are derived from 
the assumption that the sdB stars in these close bina- 
ries should rotate synchronously with the orbit. That 
means that a combination of the radial velocity, rota- 
tional velocity of the sdB star and an estimate for the 
sdB star mass, solves the system and thus gives the 
mass of the companion. The worrying aspect is that all 
systems found have low inclinations. If the companions 
are indeed neutron stars or black holes, as is implied by 
their inferred mass, their fraction is astonishingly large, 
with several objec ts in a sample of 31 (Gei e r et al . this 
volume). Indeed lYungelson and Tutukovl (j2005() pre- 
dict small fractions, and in that study neutron stars and 
black holes do not receive an asymmetric kick, which 
typically drastically reduces the surviving close binary 
population. 

Using the same population synthesis model as de- 
scribed before, including a symmetric kicks at the for- 
mation of neutron stars (see Portegies Zwart and Yungelsoii^ 



1999t Nelemans. Yungelson. and Portegies ZwartlboOl ) 
we find birth rates for sdB stars with neutron stars and 
black holes of 1.6 x IQ-'^ and 2.2 x lO^V"^- This 
implies about 1 sdB - neutron star binary per 100 sdB 
binaries and 1 sdB - black hole binary per 10,000 sdB 
binaries. It confirms that finding so many massive com- 
panions indeed is not expected, so very exciting if con- 
firmed by higher inclination systems. 



6 Conclusions 

sdB stars are a very useful population for constraining 
binary evolution theory via population synthesis stud- 
ies, because they probe very particular evolutionary 
channels and the number of known objects and their 
properties is growing rapidly. I have presented prelimi- 
nary results of a population synthesis study of sdB stars 
and related objects. Our models show differences with 



earlie r studies, s uch a s lYungelson and Tutukovl fe.g 



2OO5I): iHan et al\ (e.g. boOSl ): iTutukov and Yungelson 



(e.g. 119901 ) due to differences in assumptions about bi- 
nary evolution and the star-formation history of the 
Galaxy. This means that when proper modelling of 
the selection effects that govern the detectability of the 
observed sample is done, we can use the sdB stars to 
constrain the binary evolution models. 

I have shown the promise of comparing the popu- 
lation of sdB stars with their direct cousins, such as 
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helium-core white dwarfs and have compared the birth 
rates of (apparently) single sdB stars and (apparently) 
single helium white dwarfs. I conclude that the birth 
rate of single helium white dwarfs is substantially larger 
than that of single sdB stars, suggesting that at hast for 
the helium white dwarfs interactions on the red giant 
branch with sub-stellar companions are important. 

Finally we have calculated the birth rate of sdB stars 
with neutron star and black hole companions and con- 
cluded that the models predict substantially smaller 
fractions (1 and 0.01 per cent respectively) than found 
in short period sdB binaries, assuming the sdB star is 
in co-rotation with the orbit. This either means that 
that assumption is not correct, or much more interest- 
ingly, that there are many mode sdB stars with massive 
compact companions in close binaries than the models 
predict. 
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